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Abstract: This paper provides a class of feedback controllers that guarantee global stability of quantum 
angular momentum systems. The systems are in general finite dimensions and the stability is around an 
assigned eigenstate of observables with a specific form. It is realized by employing the control law which was 
proposed by Mirrahimi & van Handel. The class of stabilizing controllers is parameterized by a switching 
parameter and we show that the parameter between and 1 /N assures the stability, where TV is the dimension 
of the quantum systems. 



1 Introduction 

Recently it has been suggested that the stabilization of 
quantum angular momentum systems is required to re- 
alize quantum information technologies |2J. Feedback 
control laws proposed before Mirrahimi & van Han- 
del are applicable only to low dimensional quantum 
angular momentum systems [31 1131 ITS] . However, the 
stabilization of higher dimensional angular momentum 
systems is required for realizing these technologies [0] • 

In these situations, Mirrahimi & van Handel EI have 
developed a control law which globally stabilizes gen- 
eral finite-dimensional quantum angular momentum 
systems around an assigned eigenstate of the observ- 
ables. As far as we know, this is the first result which 
gives a method to control general finite-dimensional 
quantum systems. The purpose of this paper is to 
support the result of jHJ- 

The control law in 8 is parameterized by a switch- 
ing parameter 7. Mirrahimi & van Handel [H] shows 
the existence of the switching parameter 7, which 
guarantees the stability of the quantum systems. From 
the view point of designing controllers, the result of 
[5] does not make clear what value of 7 the designer 
should set for controllers. This is our motivation in 
this paper, where we will explicitly show a set of 7 
which guarantees the global stability. Thereby we can 
synthesis a stabilizing controller. 

The paper is organized as follows. In Section 2 we 
briefly review the quantum feedback control and intro- 
duce theorems used in this paper. Section 3 formulates 
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the stabilizing problem investigated in this paper and 
explains the result by Mirrahimi & van Handel [5]. 
Section 4 is devoted to the main result, where we pro- 
vide a class of controllers which globally stabilize the 
system and guarantee the convergence to an assigned 
eigenstate. A numerical example is given to show the 
effectiveness of the main result in Section 5. In Sec- 
tion 6, we conclude the paper. 

2 Preliminaries 

In this section, we briefly review the quantum feedback 
control and show two theorems used in Section 4. This 
section includes fundamental results about quantum 
mechanics and stochastic processes. Note that some of 
these results are explained with a specific form which 
is naturally applicable to quantum feedback control 
problems in this paper. 

Quantum state In this paper, we consider a quan- 
tum system with dimension 1 < N < 00. The quan- 
tum state of the system is denoted by p, an operator 
in a Hilbert space, which belongs to the state space 



S = { P e 



^NxN 



:p = p*, Trp = l, p>0}, (2.1) 



where p* denotes Hermitian conjugation of p and 
<C NxN is a set of all N x N complex matrices. 



Quantum measurement We consider the case of 
that observables are Hermite matrices in C NxN . We 
also assume that an observable A is not degenerated, 
that is A has mutually different eigenvalues {a^}. In 
the case of orthogonal measurement [7] with the ob- 
servable A of a quantum state p, the numerical out- 
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come is randomly selected from {a,k}(k — 1, . . . , N). 
The outcome a k is observed with probability 

Prob(a & ) := rp(a k )*p ^(a fc ), (2.2) 

where ip(ak) denotes a corresponding normalized 
eigenvector of A. This measurement causes a jump 
of the state p to the eigenstate of A, i/'( a fc)V'( a fc) ■ 
Because of this jump it is difficult to realize a quan- 
tum feedback control by an orthogonal measurement. 
This is a motivation for employing a continuous mea- 
surement [Tim]. 

Feedback control of quantum systems In this 
paper, we consider quantum controlled systems via a 
continuous measurement. The dynamics of such a sys- 
tem with feedback obeys the following Ito equation for 
the conditional density pt , which is a quantum analog 
of the Kushner-Stratonovich equation of nonlinear fil- 
tering PE1ITT]: 



dp t 



- i[H, p t ]dt - iu t [G,p t ]dt 

+ (cp t c* ~ ^(c*cp t + p t c*c))dt 

+ ^j(cp t + Pt c* - Tr[(c + c*)p t ] Pt )dW t , (2.3) 



where 



Pt : a quantum state at time t, 
po : an initial state, 
dWt : an infinitesimal Wiener increment satisfying 
the ltd rules : V[(dW t ) 2 ] = dt,E[dW t ] = 0, 
H : system's intrinsic Hamiltonian, 
G : a control Hamiltonian, 
u t : control input (u t G M), 
c : an observable 

1) : the detector efficiency (0 < 1) < 1). 

The third and fourth terms in (|2.3|) represent the de- 
terministic and stochastic back-action of the measure- 
ment, respectively. It should be noted that the solu- 
tion of (|2.3() is continuous in time |10|. 

We will investigate the control problem of feedback 
stabilization of the equation. 

Stability of the quantum system The solution 
of the equation l|2.3|l is a Markov process. We define 
the stochastic stability of quantum systems as below 
by following Kushner jo]: 



The equilibrium p c is said to be stable in probabil- 
ity if 

Ve > 0, 3r(e) : (0,oo) -> R+ := [0,oo) 
s.t. 

\\po - Pc|| < r(e) 



^Pr sup Upt-Pell > e =0, (2.4) 

\0<t<oo / 

where \\ ■ \\ is an arbitrary norm of a matrix in 

£NxN 

2. The equilibrium p c is globally stable if it is stable 
in probability and additionally 

V Po eS Pr f lim p t = p e ) = 1. (2.5) 

Preliminary results Before closing this section, we 
introduce two theorems which will be used in the 
proofs in Section 4. 

We can apply a theorem in 0] pp. Ill, Lemma 4.3 
to (|2.3|l from the viewpoint that Ij2.3|l is a transition 
equation of Markov process. Thereby we obtain the 
following theorem: 



Theorem 2.2. [< 

starting from pq. 



Consider a diffusion process pt G S 
Let r be an open subset of S and 



j (r) be the first exit time of pt from T. Then for all 



T>0, Pa G S, 



i 

E[Tp ° (r)] " l-sup poeS Pr{r Po (r)>T} 



(2.6) 



The other theorem is applicable to state equations 
for deterministic systems which is investigated in Sec- 
tion 4. We quote the theorem from [Sj. 

Definition 2.3 (Invariant set [H]). An invariant set C 
is defined as a set with the property that if the initial 
state of the system is in C then its whole path (forward 
and backward) lies in C . 

Theorem 2.4 (LaSalle's invariance principle 6 ). Let 
Q(p) : S — ► be a scalar function with continuous 
first partial derivatives. Let designate the region 
where Q{p) < I. Assume that tti is bounded and that 
within ill, Q(p) is positive definite and ^Jf^ < 0. Let 
R be the set of all points within i7/ where ^Jf = 0, 
and let C be the largest invariant set in R. Then every 
solution p t in fli tends to C as t — * oo 



3 Control law 



Definition 2.1. 
1(231, Le - dp t \ Pt = 



Let p e be an equilibrium point of The purpose of this section is to introduce the control 
, o = 0. law given by Mirrahimi & van Handel |H]. 
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Stabilizing control problem At first we explain 
the control problem of the quantum systems 1)2.3)1 as 
follows. The quantum systems are assumed to be un- 
der continuous measurement of a fixed angular mo- 
mentum J ( J is a positive integer or half- integer ). 
This means the dimension N of the systems is 2 J + 1. 
A typical case of such systems is an atomic ensemble 
detected through a dispersive optical probe ^21- We 
measure the angular momentum of the systems along 
an axis, say z and control the angular momentum by 
applying a magnetic field to the system along an axis, 
say y which is in a direction perpendicular to an axis 
z. The control objective is to move the quantum state 
into an eigenstate of the angular momentum. 

The dynamics of such quantum systems under an 
ideal condition can be described by 1)2.3)1 with the 
three operators; H , G, c, defined as follows: 

1. H = 0. 

2. G = [3F y (0 > 0), where F y represents the angular 
momentum along the axis y of the form [7] 



Fv = — 

y 2i 



-ci 
ci 



O 







-c 2 



CN-2 -CJV-1 

cjv-i 



C fe = ^/{N-k)k = y/{2J+l - k)k. 



(3.1) 



3. c = cxF z (a > 0), where F z represents the angular 
momentum along the axis z of the form |7| 



F z = 



-J 



O 



o 



-J + l 



J - 1 



(3.2) 



Let Afc := k — J — 1 (k = 1, . . . , N) denotes an eigen- 
value of F z and ipk is the corresponding eigenvector. 
Furthermore without loss of generality we assume a = 
P = 1 by scaling of the time and u t and we obtain: 

dp t = - iu t [Fy,p t }dt - ~[F Z , [F z ,p t ]]dt 

+ ^j{F zPt + p t F z - 2Tr[F z p t }p t )dW t7 (3.3) 

where < r) < 1. The control objective is to stabi- 
lize ,0(y) := V'/V'/ globally where / is one of indices 
1, . . . , N . We define a distance function 

V(p) := 1 - Tt( PPu) ) = 1 - (p) ff :S^[0, 1] (3.4) 

from the state p to the target state for the expla- 
nation of the result of [8] and the proof of our main 
result. 



Mirrahimi & van Handel 8 have proposed a control 
law for this stabilization problem as in the following 
theorem: 

Theorem 3.1. jS] Consider the system 1)3. 3|l evolving 
in the set S. Let 7 > 0. Consider the following control 
law: 

1. u t = -Tr(i[F y ,pt\p(f)) if V(p t ) < I-7; 

2. u t = 1 if V(p t ) > 1 - 7 /2; 

3. If Pt e B = {p : 1 - 7 < V(p) < 1 - 7/2}, 
then u t — —Tr(i[Fyip t ]p(f)) if Pt last entered B 
through the boundary V(pt) = 1 — 7, and Ut = 1 
otherwise. 

Then ^7 > s.t. u t globally stabilizes ()3.3(l around 
ptf) and E[p t ] — > p^^ as t — > 00. 

Hereafter we call the control law in Theorem 13.11 
MR control law and «mh(t) denotes the control input 
using MH control law with the parameter 7. 

4 Synthesis of stabilizing con- 
trollers 

4.1 Main theorem 

Theorem 13 . II shows the existence of a positive parame- 
ter 7 with which the stability of the quantum systems 
is guaranteed. In other words, Theorem 13. II is the ex- 
istence theorem for the stabilizing controllers. On the 
other hand, from the viewpoint of controller design, it 
is unclear what positive number 7, which the designer 
should set for the controllers, guarantees the stabil- 
ity. In this section, we give a set of 7, which definitely 
guarantees the stability of the quantum systems, in 
the following theorem which is a main result of this 
paper. 

Theorem 4.1. Consider the system (|3.3|) evolving in 
the set S. Let u t = umh(7)- If 7 G (0,^), Ut glob- 
ally stabilizes 1)3. 3|) around the target state pit\ and 
E[pt] — > P(j) as t — > 00. 

The proof of Theorem l4.1l is executed in the follow- 
ing three steps: 

1. At first, we show that the ensemble average of the 
solution of 1)3.3)1 with a p t -independent nonzero con- 
trol input Ut € C 1 converges to the maximally mixed 
state as t — > 00. 

2. Secondly, we prove the following lemma, the proof 
of which is given in subsection 4.2: 
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Lemma 4.2. Consider the system 1(3.3(1 with Ut = 1 
and an initial state po G iS>i_ 7a := {p G S : 7 a < 
< 1}- iei T Po (5>i_ 7a ) 6e the first exit time of p 
from 5>i_ 7a . Then, i/0 < 7 a < 1/iV, 

sup £[T P0 (S>i_ 7a )] < oo. (4.1) 

Po£5>i — , a 

Remark 4.3. This lemma gives a set of 7 a which 
satisfies (|4.1() . while Mirrahimi & van Handel 8 shows 
the existence of such 7 a . This is a main defference 
between [H] and this paper. 

3. Finally, we complete the proof of Theorem 14.11 by 
combining the results of the second step and the fol- 
lowing lemma which we collect lemmas in [H] into: 

Lemma 4.4. [S] Consider the system (|3.3|l with 
umh(7) and an initial state po G S. if 7 = 7 a satisfies 
(|4.1Jl . «mh(7) globally stabilizes (|3.3(l around pin and 
E[pt] — ► /?(/) as i — * 00. 

4.2 Proof of Theorem HTTl 

This subsection is devoted to the proof of Theorem l4.ll 
which is devided into three parts corresponding to the 
three steps proposed in the previous subsection. 

Step 1. In this step, we consider the dynamics of 
E[pt] =: pt, 

dpt = -iut[F y ,pt]dt- i[F„ [F z ,p t ]]dt, (4.2) 

where a control input is nonzero and /^-independent 
and has continuous first derivatives. We call 1(4.2(1 the 
ensemble dynamics. 

At first we give the following lemma on the equilib- 
rium points of 14.2(1 . 

Lemma 4.5. The ensemble dynamics 1(4.2(1 has a 

unique equilibrium point j?I> which is the maximally 
mixed state of (|2.3(l . 

Proof. Let us assume that the ensemble dynamics has 
an equilibrium point p c . That is 

~iu t [F y ,p c ]dt~ l[F„ [F z ,p c ]]dt = 0. (4.3) 

Multiplying both sides of (|4.3|l by 2p c and subse- 
quently taking trace of the both sides gives 

Tr([F z ,[F Z7 ff]}(f)dt^0. (4.4) 



Remember F z and p e are Hermitian matrices. The left 
side of (|4.4|) can be computed as 

Tr([F z ,[F z ,p c W)dt 

= Tr{(f?yF;iF z ,p c ] - {F z ,f?}F;(f?y}dt 

= Tr{(f?yF;iF z ,p c } - F;(pr[F z ,p c }}dt 
= Tr([F Zl pr[F z ,p c ])dt 

= \\[F z ,fF}\\ 2 dt, (4.5) 

where Hfi^jp ]!! denotes the Frobenius norm of the 
matrix [F z ,p c ]. I|4.3|l and l|4.5|l show 

([F z ,p c }) lm = 

(I,m6N, 1 < l,m < N), (4.6) 

where ([F z , p°])im stands for the (Z,m)th element of 
the matrix [F z , p c ] . Noting that F z is a diagonal ma- 
trix, we obtain 

= (p c ) lm {(F z ) u - (F z ) mm } - 0, 

(l,m GN, 1 < l,m < N). (4.7) 

Note also that F z has no repeated diagonal entries. 
By ()4.7|l . p c should be a diagonal matrix. 

Now we substitute a diagonal matrix p t — p° into 
(|4.2|) and obtain 

(-i[F y ,p c ])i m 

= -i{F y )i m {{p c )mm - {P c )u} = 0, 

(l,m G N, 1 < l,m < N). (4.8) 

Since the location of the nonzero elements of F y is 
limited as in l|3.1(l and with l|4.8() , the all diagonal ele- 
ments of p c are known to be equal. Such a matrix in S 
is only the maximally mixed state jjl. This completes 
the proof. □ 

We can derive the following proposition from this 
lemma. 

Proposition 4.6. Consider a p t -independent nonzero 
control input u t G C 1 . For any initial state po G S, the 
solution pt of H4.2fl converges to the maximally mixed 
state jfl as t — > 00. 

Proof. Consider a function 

Q(p t ) := Tr(p 2 t ) - 1 (4.9) 

as a candidate of the Lyapunov function. It is easily 
verified with Tip t = 1 that Q(pt) > for all pt G S and 
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that Q(pt) = iff pt = jfl- A computation similar to Therefore, in the case e = jt — 7a, we obtain 



that shown in 1(4. 3 j) and 14.5J1 gives 
dQ(pt) 



dt 



\[f»(hW<o. 



(4.10) 



We have equality in l|4.10|l iff p t is a diagonal matrix 
(from the same argument in Lemma l4.5|l . Let C be 
the invariant set contained in {p S S : d ®Jt — 0}. We 
can see the set C includes any diagonal matrix. 
Recall the result and proof of Lemmma 14.51 The 



3Ti < oo s.t. t > Ti => E[V(p t )] < 1 - 7a- 

With this lemma, we can prove Lemma f4. 21 
Proof of Lemma\4-S\ By Lemma 14.71 



(4.15) 
□ 



maximally mixed state satisfies 



dQ(P) 
dt 



and 



it is an equilibrium point of the ensemble dynamics. 
So the maximally mixed state is an elemnt of C 

On the other hand, let us assume that p is a diagonal 
matrix other than the maximally mixed state at time 
t. We can obtain 

(dp t )lm = -iUt(Fy)lm{(Pt)mm ~ (p~t)ll}, 

(l,m e N, 1 < l,m < N,l ^ m, ). (4.11) 



In the equation (|4.11() . there exist I and m such 
that {dptjim is not zero as we saw in the proof of 
Lemma 14.51 That is, p changes into an off-diagonal Step 3 
matrix. Therefore a diagonal matrix other than the 
maximally mixed state is not an element of C. 

Thus the maximally mixed state is the unique ele- 



3T ; T = inf {t \ E[V(p t )} < 1 - 7 a}- (4.16) 

t<E[0, oo) 

In addition, applying Theorem 12.21 to the case T = 
<S>i_ 7a and T = Tq yields to 

T 

Elr^S^J] < 1 _ 8upces p r{Tc ° (5>1 _ 7j>T()} - 

(4.17) 

Now we can prove this lemma in the same way with 
()4.16|) . (|4.17|) as the proof of Lemma 4.6 in Mir- 
rahimi & van Handel [Sj. □ 



ment of C. The assertion that pt —> jjl as t — > oo 
can be proved by applying Theorem= 12.41 (LaSalle's 
invariance principle). □ 

Step 2. In this step, we consider the system (13.311 
with u t — 1 and an initial state po G S. At first, we 
show the following lemma. 

Lemma 4.7. For 7 a £ (0, -k), there exists a finite 
time T\ such that 

min E[V(p t )} < 1 - 7a . (4.12) 

te[o, Ti] 

Proof. The solution of (|3.3|) is continuous in t and it 
is obvious that V(pt) is continuous in p t . Thus Propo- 
sition 14.61 implies 



lim E[V(p t )] = V lim E[p t ] 



(4.13) 



That is, 

Ve>0, 3Ti 

s.t. t > T, => \E[V( Pt )} - (1 - 1)| < 6. (4.14) 



Proof of Theorem \4-l\ Consider the system H3.3(l with 
the control input umh(7) and an initial state po G 
S. When po e 5>i_ 7 , the control input is Ut = 1 

until p t exits from 5>i r Therefore we can show 

that 14.1(1 holds for 7 = 7 a £ (0,^) using Lemma 
14.21 Then, applying Lemma [4.41 completes the proof 
of Theorem O □ 



5 Numerical example 

We here give a numerical example to show the effec- 
tiveness of the main result in this paper for synthe- 
sis of stabilizing controllers. We deal with a quantum 
system under a continuous measurement of a fixed an- 
gular momentum J — 10, that is the dimension of the 
system N is 21. Let us stabilize (|3.3ll around pni) 
globally from the initial state p^ by utilizing the MH 
control law. 

We can see from Theorem 14. II that we only have to 
set a switching parameter 7 between to = ^ in 
order to realize the global stability. Figs. [J and [5] il- 
lustrate three sample paths of V(p t ) in the case that 
7 is set to 0.04(< Jj-, Fig. [I} and 0.4(> Fig. EJ) 
Respectively. The value of 7 = 0.04 guarantees the 
stability of the 21-dimensional quantum angular mo- 
mentum system. However, the value of 7 = 0.4 does 
not so. These numerical simulations illustrate the re- 
sult of Theorem 14.11 and show the effectiveness of it. 
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Figure 1: In the case of 7 = 0.04 < all the sample 
paths converge to pnx) from pn)- 




Figure 2: In the case of 7 = 0.4 > A there exists a 
path which is expected not to cnverge to P(u)- 

6 Conclusion 

This paper provides a class of controllers that guar- 
antee the global stability of an assigned eigenstate for 
./V-dimensional quantum angular momentum systems 
by utilizing the result of [H] ■ The main result of this 
paper is in Theorem 14. II which says that any switch- 
ing parameter 7 G (0, A) in a MH control law given 
by [5] is sufficient for stability in probability 1. 

However, we should note that Theorem 14.11 gives 
only a sufficient condition for the stability in proba- 
bility 1. A necessary and sufficient condition for the 
stability remains as a future interesting work. 
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